arXiv:1508.07284v2 [quant-ph] 20 May 2016 


Loschmidt echo in many-spin systems: contrasting time-scales of local and global 

measurements. 

Pablo R. Zangara,^ Denise Bendersky/ Patricia R. Levstein/ and Horacio M. Pastawski^ 

^Instituto de Fisica Enrique Gaviola (CONICET-UNC) and Facultad de Matemdtica, 

Astronomia y Fisica, Universidad Nacional de Cordoba, 5000, Cordoba, Argentina 

A local excitation in a quantum many-spin system evolves deterministically. A time-reversal 
procedure, involving the inversion of the signs of every energy and interaction, shonld produce the 
excitation revival. This idea, experimentally coined in NMR, embodies the concept of the Loschmidt 
echo (LE). While such an implementation involves a single spin autocorrelation Mip, i.e. a local LE, 
theoretical efforts have focused on the study of the recovery probability of a complete many-body 
state, referred here as global or many-body LE Mmb- Here, we analyze the relation between these 
magnitudes, in what concerns to their characteristic time scales and their dependence on the number 
of spins N. We show that the global LE can be understood, to some extent, as the simultaneous 
occurrence of N independent local LEs, i.e. Mmb ~ This extensive hypothesis is exact 

for very short times and confirmed numerically beyond such a regime. Furthermore, we discuss 
a general picture of the decay of Mi,i as a consequence of the interplay between the time scale 
that characterizes the reversible interactions (T2) and that of the perturbation (te). Our analysis 
suggests that the short time decay, characterized by the time scale te, is greatly enhanced by the 
complex processes that occur beyond T2 . This would ultimately lead to the experimentally observed 
T3, which was found to be roughly independent of te but closely tied to T2. 
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I. INTRODUCTION 

If an ink drop falls into a pond, the stain diffuses away until no trace of it remains whatsoever. One may naturally 
say that such a process is in fact irreversible. In the microscopic world, similar phenomena are also ubiquitous. For 
instance, let us consider a many-spin quantum system in thermal equilibrium where a local polarization excess is 
injected. Then, this excitation would spread all over as consequence of spin-spin interactions. Such an apparently 
irreversible process is known as spin diffusion 0,i and it can lead the system back to equilibrium. However, this 
naive picture has its limitations. On the one hand, spreading is not always the rule as there are physical situations 
where the initial excitation does not vanish. This is the case of Anderson localization [^, or when the excitation 
remains topologically protected On the other hand, even in the cases where the system seems to have reached 
an equilibrium state, the unitarity of quantum dynamics ensures a precise memory of the non-equilibrium initial 
condition. Then, if some experimental protocol could reverse the many-body dynamics, it would drive the system 
back to the initial non-equilibrium state Q. Such a general idea defines the Loschmidt echo (LE), which embodies 
the various time-reversal procedures implemented in nuclear magnetic resonance INMR) 

The hrst NMR time-reversal experiment was performed by Hahn in the 1950’s [l0|. The procedure, known as spin 
echo, reverses the precession dynamics of each independent spin around its local magnetic field by inverting the sign 
of the Zeeman energy. However, the sign of the energy associated to the spin-spin interactions is not inverted and, 
accordingly, the echo signal is degraded. Such a decay occurs within the time scale T 2 that characterizes the spin-spin 
interactions. Indeed, these interactions determine the survival of a spin excitation at short times as ^ 1 — and 

its later complex dynamics generating a diffusive spreading. By the early 1970’s, Kessemeier, Rhim, Pines and Waugh 
implemented the reversal of the dynamics induced by the spin-spin dipolar interaction [111 Il 2 | . This results in the 
“magic echo” which indicates the recovery of a global polarization state. Two decades later, Ernst and collaborators 
introduced the “polarization echo” [l3l| . There, a local excitation injected in a many-spin system is let to evolve, 
then time-reversed and hnally detected locally at the initial spot. While the success of these time reversal echoes 
unambiguously evidenced the deterministic nature of spin-dynamics in NMR, it is clear that the reversal is unavoidably 
degraded by uncontrolled internal or environmental degrees of freedom or by imperfections in the pulse sequences. 
Furthermore, the degradation seems to occur in a time scale, say Ta, much shorter than a naive estimation of the 
characteristic scale of these perturbations, say ts. Then, the question that arises is whether the complexity inherent 
to a large number of correlated spins would enhance the fragility of the procedure under perturbations. 

A next generation of experiments in organic crystals |14l - ll6l | seemed to confirm that the experimental never 
exceeds a few times T 2 . In other words, T 3 keeps tied to the time scale that characterizes the reversible many- 
body interaction. This led to postulate that in an infinite many-spin system the complex dynamics could favor 
the action of any small non-inverted interaction that perturbs the reversal procedure. Thus, reversible interactions 
become determinant for the irreversibility rate. This constitutes our Central Hypothesis of Irreversibility. Such a 
wisdom is further reinforced by the natural association of many-body complexity with a form of chaos [l^ . [l8l| and 
the conhrmation that quantum dynamics of classically chaotic systems should manifest a dynamical instability [l9j | 
which leads to an environment-independent decoherence rate [^. . 

During the last decade, solid-state NMR kept providing a versatile testing bench to study time reversal in large spin 
arrays [22l426l| . As a matter of fact, a standard experiment involves a crystalline sample with an infinitely large number 
of spins. In contrast, the numerical test of many spin dynamics has to be restricted to strictly finite systems [13, 
While this appears to be a major limitation, it allows the analysis of a situation that the experiments cannot achieve: 
moving progressively from small systems to larger ones with a fully controlled perturbation. The expectancy is that a 
sort of finite size scaling may allow to identify an emergent mechanism that rules reversibility in the thermodynamic 
limit [ 2 ^. As in the experiments, the witness for such a transition should be the LE as measured by a single spin 
autocorrelation function Mi.i, i.e. a local polarization. For short we call Mi 1 the local LE. It is not difficult to probe 
that Hi^i = (Mi l + l )/2 is the probability that a given spin, say the 1 ^*, remains up after the whole procedure. 
Besides, in a case of a spin excitation in a ID chain with XY interactions jsol - fs^l Mi 1 precisely coincides [131 with 
the global overlap of two one-body wave functions as defined for semiclassical systems [13, IH, . The square of 
the overlap between the initial and final many-body wave functions, Mmb, dehnes a global or many-body LE. It is 
important to notice that Mmb has not been addressed experimentally, but nevertheless it is a natural magnitude in 
theory (33 - [33| . Thus, we are left without a precise relation between the object of theoretical studies and experimental 
ones, i.e. Mmb and Mi^i respectively. This missing link is the central question we address in this paper. 

Here, we consider a system of N spins whose initial state is given by a local excitation injected in an high temperature 
state. Firstly, we discuss the formal relation between Mmb and Mi^i, which is derived exactly at least for very short 
times. In particular, we assess how the A^-dependence or extensivity of Mmb is evidenced in the time scales involved. 
This leads us to hint that the revival of a many-spin state results from the recovery of each individual spin configuration, 
much as if they were statistically independent events. Since in the initial high temperature state there are N/2 spins 
up, their rough statistical independence would lead to a behavior of the sort of Mmb ^ ^ (Mi_i)^^^. This 
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is confirmed by the numerical evaluation of the LE in a specific spin model. 

Furthermore, we discuss a general picture beyond the short-time regime, where the decay of Mi i results from the 
interplay between the time scale that characterizes the reversible interactions (T 2 ) and that of the perturbation (rs). 
This would ultimately lead to the experimentally observed T 3 . In such a sense, our analysis provides a conceptual 
hinge between the theoretical and the experimental realms. 

The paper is organized as follows. In Sec. Ullwe introduce the LE framework: the initial state and the time-reversal 
procedure. Here, we define both the local and global LE. In Sec. IIIII we compute the short time expansions for the 
local LE and its non-local contributions (in particular, the many-body LE). This allows us to discuss a general picture 
of the LE decay in terms of the times scales that characterize Mmb and Mi.i. The dependence with N is discussed 
in terms of the extensivity of Mmb and statistical independence of the local autocorrelations. In Sec. IIVI we assess 
our expectancies by a numerical evaluation of the LE in a spin system. Section |V] summarizes our main conclusions 
and some of the important open questions in the field. 


II. THE LOSCHMIDT ECHO IN SPIN SYSTEMS. 

Let us first specify the initial condition of a “local excitation in many-spin system”. We consider N spins 1 /2 in 
an infinite temperature state, i.e. completely depolarized mixture, plus a locally injected polarization, 

p^ = ^(X + 2Sl). ( 1 ) 

Here, the spin 1 is polarized while the others are not, i.e. tr\S2^Po] = Such an initial state can be experimentally 

implemented not only in NMR but also in cold atoms [s^. 

As in the early LE experiments HQ, our numerical evaluation focuses on an imperfect time reversed evolution of 
the excitation, followed by a local measurement. The procedure is depicted in Fig. [TJ A many-spin Hamiltonian Hq 
rules the forward evolution of the system up to a certain time t/j. At that moment, an inversion of the sign of Hq is 
performed, leading to a symmetric backward evolution. Nevertheless, there are unavoidable perturbations, denoted by 
E, that could arise from the incomplete control of the Hamiltonian, acting on both periods. Thus, evolution operators 
for these t^-periods are = exp[—and U_{tR) = exp[—^(—IL q + respectively. It is quite 

practical to define the LE operator as: 


lJLE{‘^tR) = lJ_{tn)U+{tR), ( 2 ) 

which produces an imperfect refocusing at time 2tR. A local measurement of the polarization, performed at site 1, 
defines the local LE: 


Mia(t) = 2tr[SlU^E{t)pQUl^{t)\ = 2tr[Slpt]. (3) 

Here, we choose as free variable t = 2tii, the total elapsed time in the presence of the perturbation. The time 
dependence of pt in the Schrbdinger picture is, 

Pt = Ui^E{t)pQUlj,{t). (4) 

Using equation o, and after some algebraic manipulation, the LE can be explicitly written as a correlation function: 


MiAt) = ^tr[ulAt)smuLEA)km 


tr[smsm] 

tr[^f(0)5f(0)]' 


Here, the time dependence is written according to the Heisenberg picture. 


( 5 ) 


= uiE{t)smuEEit). ( 6 ) 

Notice that equation ([5]) is an explicit correlation function at the same site but different times, i.e. an awtocorrelation. 
This kind of a correlation has been recently employed to address localization phenomena in spin systems [1, Q, |4l| 
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FIG. 1 . The scheme of the local LE, (a) The initial state is given by a local excitation in a high temperature spin 

system, as stated in equation ©. The system is let to evolve under the Hamiltonian Hq + E and the excitation diffuses all 
over until a time t = ta (b). At that time, a time reversal is performed, leading to a backward evolution ruled by —Hq + E. 
At time t = 2tR (c), a local measurement is performed at the initial spot. 


and it generalizes the standard one employed to assess spin diffusion [d^. In terms of the Hilbert-Schmidt inner 
product between the initial and the time evolved density matrices, i.e. equations o and dH) respectively, the LE can 
be written as [3,13 S] : 




2^tr[pQpt\ -1 = 2 


tr[poPt] 

tr[poPo] 


- 1 , 


( 7 ) 


which, in the present case, progressively decays from 1 to 0 as it occurs with the statistical overlap between two wave 
packets in the standard LE definition (20|. 

Equivalent expressions for the LE can be derived decomposing the statistical state in a simpler basis. In order to 
proceed with the pure state decomposition of po, we consider the computational Ising basis {|/3i)}i also known as 
5'^-decoupled basis. Additionally, we define the set A of indexes j that label basis states which have the 1®* spin 
pointing up, i.e. j ^ A \ j5j) = \Pj)- It is straightforward to verify that po = Then, 

as introduced in Ref. [3| . 


=2 


1 




2 1 


1' 

~ 2 

= 2 

ni,i(t) - 2 


( 8 ) 


Here, ni.i(t) denotes the probability that the I'** spin keeps pointing up after a time t. After some manipulation. 
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Wj \ ft) 

-E (/3,it/Li5Wift) 

)l 


(9) 


Here, B stands for the complement of A, i.e. j £ B A- Sf \(3j) = —^ \(3j). One can naturally identify and define 
the two terms that contribute to the local polarization Mi.i(t). The first sum in equation ([S]) stands for the average 
probability of revival of the many-body states, denoted by Mm sit), 
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Mm sit) = 


ieA 


2N-1 




( 10 ) 


The second sum in equation @ represents the average probability of changing the configuration of any spin except the 
1®*. The third sum stands for the average probability that the 1®* spin has actually flipped, i.e. of all those processes 
that do not contribute to Then, the processes that contribute to but not to MMB{t) are denoted as: 


Mx{t)=Y, 


2N-1 




i&A \j&A (j^i) j&B 

This balance of probabilities leads to the appropriate asymptotic behavior of 
that constrain the evolution. The identification 



according to the symmetries 


= MMB{t) + Mx{t) 


( 12 ) 


is a crucial step for the following ^discussions. 

If we use the identity SI = Sf in equation ([S]), the invariance of the trace under cyclic permutations ensures 

that (t).5f (0)] = (O)S'f (t);5(^(0)] — ^tr[Sf{t)]. Since tr[Sf{t)] = tr[.Sf(0)] = 0, then: 


Mi,i(t) = 2E ^ (Ai s^mUt)smuLE{t)st{o) m 

i 

= 2E ^ UlEit)S!ULEit) m , (13) 

iGA 

which is indeed an explicit way to rewrite equation ([3]) in the form of an ensemble average. Remarkably, since S'f is 
a local (“one-body”) operator, its evaluation in equation m can be replaced by the expectation value in a single 
superposition state (i^ . 


Mia(t) = 2 ('k„e,| UlE{t)S!ULEit) I '^neq) ? (14) 

where: 

IW = E^^e^^^ m. (15) 

Here, is a random phase uniformly distributed in [0, 27r). As a matter of fact, the state defined in equation (ITSl) is 
a random superposition that can successfully mimic the dynamics of ensemble calculations and provides a quadratic 
speedup of computational efforts [45l - l47j . 


III. THE LOSCHMIDT ECHO DYNAMICS 
A. Short time expansions and beyond. 

In order to analyze the A-dependence of the LE and its time scales, we compute here the short time expansion of 
the magnitudes MMB{t) and Mx(t). Up to 2"^^ order in time, 

Mia(t = 2tR) = 2 E ^ m UlEit)StUEEit) |ft) 

iGA 

= 1 - [t/Kf E ^ ((ft I Ift) - 2 (ft I I A)) + o (it/hf) . 

iGA 


( 16 ) 
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Similarly, the leading contributions to MMS^t) and Mx{t) are: 




i&A 
= 1 - 


(VS)" ^ ^ ((A IS" W^) -im I A)") + O [(t/hf) 

iGA 


(17) 


and 


Mx{t) = E ^ E ^LEit) IA)f - E |(/3.l U^Eit) lA) 

ieA \jeA U^i) les 

= (VS)" E ^ (2 (Al m - (ftl s lA)") + o {{t/hf). 

i&A 


(18) 


Let us consider a generic secular (i.e. polarization conserving) perturbation E given by a Hamiltonian with 
arbitrary anisotropy a, 


an 


E = 


N 


(19) 


This is still quite general since even a double quantum perturbation {S^S'^ + A A ’ 'which does not conserve 
polarization) can be reduced to a secular one by the truncating effects of radiofrequency fields [ 2 ^. In addition 


do not consider here the case 


E,A 


= 0 (e.g. pure Ising or on-site diagonal disorder), since in such a condition the 


first non-trivial order in time is the 4*^ (see Appendix). Then, the following identities hold: 

E 2 JV -1 (^*1 ^ 2Na‘^ (l~ g) ’ 

i&A ^ '' 

E jA I-®'’ = (v + 


2 ^ 


i&A 


( 20 ) 

( 21 ) 

( 22 ) 


Here, stands for the average local second moment of E, 


N 


N 


TV E AT E 


2=1 




2=1 


iV 






V 2 


(23) 


Being a perturbation, we notice that cr^ is much smaller than the average local second moment ctq of the unperturbed 
Hamiltonian Hq. In terms of time scales. 


A = fi/y^ < h/ = te. 
The identities in equations (l20ll . (1211) and (l2^ lead to 




(24) 


(25) 


and, 


MMB{t) = 1 - JiV (t/rs)" + O (it/hf) , 

MT,(t) = {t/r^f + O (it/hf) . 


(26) 

(27) 
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These expansions hold for t < (ts/N). Beyond such a very short time regime, a general term in the expansion of 
Ml lit) will be of the form 


CiN.n)t^/ {rhT^-^) (28) 

with k > 2 and the coefficient C( 7 v,n\described by combinatorial numbers of increasing size that depend on the 
topology of the interactions (e.g. see 13)- Since the experimental set up corresponds to the limit described by 
equation (I24L this expansion will be dominated by terms with the lowest possible order in the weak interaction, i.e. 
k = 2-. 




n 


(29) 


Equation (|2^ indicates that beyond the very short time expansion, i.e. (r^/iV) < t < the dependence on 
becomes superseded by the diverging terms in the scale T 2 . This could lead to the new time scale T 3 which was seen 
experimentally to be tied to T 2 as 


T 2 <T 3 <Te. (30) 

In that sense, T 3 becomes characteristic of the complexity or “chaos” of the many spin system that amplifies the 
small effect of the perturbation. In addition, it is important to stress that, being an experimental fact, equation (1301) 
corresponds to a system composed by infinitely many interacting spins. In other words, equation (1301) stands for the 
relations of time scales in the thermodynamic limit. Quite on the contrary, any numerical simulation involves a finite, 
very small indeed, number of spins where the irreversibility rate T 3 would be essentially given by ts- Then, the LE 
decay rate evaluated in a finite system would ultimately be perturbation-dependent [ 2 ^. Thus, our Central Hypothesis 
of Irreversibility would mean that equation (1301) is an emergent property. It should rely on the thermodynamic limit, 
which implies taking the limit N ^ 00 first, and then —>■ 00 . The non-uniformity of these limits plays a crucial 
role to yield quantum phase transitions, as discussed in the context of Anderson localization [49l - [^ . 

The physical picture described above is schematically represented in Fig. [51 There, we show the expected interplay 
between Mm sit) and Mx{t) leading to Mi^i(t). Indeed, as stated in equations ([55)) and (1571) . the very short time 
dependence of both contributions is extensive in N: Mmb(Q decreases as 1 — Na'^t^/A and Mxit) increases as 
[N — 4)CT^t^/4. Such a precise balance provides for the short time decay of Mi^iit) given by equation (1551) . i.e. 
1 — Notice that there is no reason to assume that the decay of Mm sit) would remain ruled by r^. Beyond 

the very short times, we expect that the time scale T 3 should also show up as T^/N'^ with ^ 1 in the decay of 
MMBit)isee below). Furthermore, while Mm sit) goes monotonically to zero, Mxit) displays a highly non-trivial 
behavior. Indeed, Mxit) first increases by feeding from the decay of Mm sit) until it reaches a maximum. This 
growth indicates a progressive divergence of long-range correlations. Thereafter, Mxit) should decay accounting for 
the fact that the state remains properly normalized. This is precisely what Mi^iif) measures: a conserved polarization 
that ultimately distributes uniformly within the spin system. In an isolated finite system this implies the asymptotic 
plateau Moo 1/N. As pointed above, the decay of both Mi^i(t) and Mxit) occurs in a time scale T 3 , which 
according to equation (1551) . is somewhat longer but close to the “diffusion” time T 2 . This is the regime captured 
experimentally. 


B. The extensive decay hypothesis 

The previous short-time expansions provide a hint on the scaling relation between the local LE, Mi_i(t), and the 
global one as embodied by Mm sit)- In particular, let us first compare the probability of refocusing the configuration 
(up or down) of a single spin, i.e. ni^i(<), and the probability of refocusing a complete many-spin state Mmb(Q- If 
the refocusing of each individual spin could be treated as an independent event, then the scaling between Hi 1 and 
Mmb would be extensive in N, 


iHi^iit))^/^ ^MMBit). (31) 

Here, the factor 1/2 in the exponent comes from equation (|T]), i.e. the initial high temperature state, where basically 
half of the spins point up, and half of them point down. Then, one can resort to the picture of a lattice gas where 





Time 


FIG. 2. A pictorial scheme of the time dependence of (dotted line) and its contributions, i.e. Mm sit) and Mx{t) (solid 

lines). Their short-time expansions, as stated in equations (1251) . (1261) and (1271) . are indicated with arrows. In particular, the 
expansion corresponding to the short-time behavior of Mi,i(t) is plotted with a dashed line. (Online version in colour.) 


N/2 particles jump among N lattice sites. As in the well known Jordan-Wigner transformation [5^, a fermion is 
associated to a spin pointing up and a vacancy corresponds to a spin pointing down. Thus, the microstate of the gas 
is completely described by the position of N/2 particles. 

Strictly speaking, the notion of extensiveness corresponds to standard thermodynamic quantities such as the entropy 
of the system. In addition, as discussed in 0 , S = —ln{Mi^i{t)) is precisely a measure of the entropy. Then, the 
validity of equation m implies an extensivity relation between the entropy per spin and the total entropy of the 
system. 

According to equations ([5]) and (1^ . 


ni,i(t) = l-i(t/rs)Vo((t/ft)'), (32) 

which in turn, up to 2““^ order in time, implies ni^i(t) ~ (Mi_i(t))^^^. Thus, equation (1^ yields 

(33) 

This is precisely the relation verified between equations (1^ and (051) . 

One might expect that beyond the very short-time decay, individual spin autocorrelations deviate from the statistical 
independence. However, this deviation will still have a local nature and therefore the A^-extensivity would remain 
valid. Indeed, we propose 

(Mi,i(t))''~MMB(i), (34) 

where the exponent rj would be some appropriate function 77 = r]{N,t). Our “extensive decay hypothesis” implies 
that 77 factorizes: 

77 (A,t)=Ax/(t), (35) 

where f{t) stands for a function that encloses information of the correlations originated by the system dynamics. 
Additionally, 


lim fit) 

t-yO+ 


1 

4' 


is required in order to recover equation (1551) . i.e. the statistical independence. 


(36) 








9 


IV. A ID MODEL. 

The physical picture described above is discussed here under the light of a specific model. In particular, we assess 
the validity of equations dSSl) and (1551) . We consider a 1-D spin chain with an anisotropic interaction described by: 


N-l 


2 = 1 


Ho = J2Jo( + S^S- 




cy qy 
*^ 2+1 


(37) 


with periodic boundary conditions, i.e. a ring configuration. Here, Jq stands for the natural units of the spin-spin 
interaction energy. As a perturbation E we choose a next nearest neighbors interaction described by: 


N-2 

i=l 

Such a perturbation appears naturally when one attempts to build an effective one-body dynamics from linear crystals 
with dipolar interactions [h^ . This is also the case in a regular crystal, when the natural non-secular dipole-dipole 
terms are truncated by the Zeeman energy of the radiofrequency irradiation, which ultimately leads to effective secular 
two-body next nearest neighbors interactions [2^. 

The local second moments and CTq of E and Hq respectively can be evaluated as in equation (|2^ : 


qZ QZ I QX QX , qy qy 

2^2 *^2 + 2 "t“ ^2 + 2 "T 


(38) 



(39) 


(40) 


and constitute the main energy scales of our problem. 

In Fig. [3] we plot MMB{t) and Mx{t) for the particular choice Je = O.IJq. Short-time expansions given 

in equations (1^ . (1^ and (l?7l) are evaluated according to equation (15^1) . It is observed that MMB{t) vanishes for 
long times. Actually, a close observation shows that MMB{t —>■ oo) ^ 0{2~^) (data not shown). In addition, notice 
that Mx{t —>■ oo) ^ 1/N. Such an asymptotic contribution provides for the equidistribution of the spin polarization 
—>■ oo) ^ 1/N. This long-time saturation corresponds to the equilibration of a finite system. 

In contrast to our schematic plot in Fig. [51 here Mx{t) does not get too close to I and MMB{t) does not decay 
much faster than Mi^i(t). Since Mx{t) provides for the whole once that MMB{t) has fully decayed, the 

contribution of Mx(t) is considerable only at long times. These effects are a consequence of the relatively small size 
of the system considered. Indeed, the case in Fig. [3] corresponds to V = 14 spins, and thus the exponent that relates 
Mi.i(t) and MMB{t) is quite small {N/A) = 3.5. The need for larger systems indicates that revealing the dominant 


orders in equation (I29|) is a major numerical challenge that may go beyond the state-of-the-art techniques 541. 


In order to assess the accuracy of the “extensive decay hypothesis”, in Fig. [4] we address the scaling relation between 
and AlMB^t) discussed in Sec. IIIIBI In particular, we try out the factorization stated in equation (1551) . By 
plotting \og{MMB{t))/ (log(Mi_i(t))V) as a function of time, we observe a unique function which does not depend 
on N 01 Je, but it has a weak dependence on time. Such a unique curve is indeed f{t) as defined in equation (1351) . 
This means that the extensivity relation between and MMB(t) is confirmed. The statistical independence, in 

turn, fails progressively once /(t) departs from the 1/4 factor of the ideal relation in equations (I55|) and (I5S1) . Since 
beyond the short-time regime f{t) decreases with time, we conclude that the recovery of a single spin is tied to the 
recovery of its neighbors. Thus, the spins are positively correlated and the revival probability of the complete V-spin 
state is enhanced. This argument is particularly relevant in ID systems. 

After the onset of the saturation regime, where ^ 1/7V and Mmb ^ the universal scaling naturally 

becomes noisy and the curves for different N and Je separate each other. Since the decay is faster for larger 
perturbations, the appearance of such a spurious behavior is observed to occur first for the largest value of Je 
considered (Je = O.SJq, plus signs and triangles). 


V. CONCLUSION 

We presented a detailed analysis of the LE in interacting spin systems. As in the NMR experiments, a local version 
of the LE, Mi^i, is defined as a single spin autocorrelation function. Simultaneously, we define a global LE, Mmb, 
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FIG. 3. The local LE and its non local contributions. MMB{t) and Mx{t) correspond to the solid lines as indicated 

by the labels in the figure. N = 14, Je = O.IJq. The short-time expansions given in equations ([251). pel) and ([27]) are shown 
in dashed lines. (Online version in colour.) 



FIG. 4. The relation log(MMB(t))/log(Mi,i(t))) as a fnnction of time. For Je = O.IJq the sizes plotted are: = 10 

(circles), A^ = 12 (squares). A'’ = 14 (diamonds), and Ai = 16 (stars). For Je = 0.2Jo the sizes plotted are: A^ = 10 (solid line), 
A'^ = 12 (dashed line). A'’ = 14 (dash-dot line), and A'^ = 16 (dotted line). For Je = O.SJq the sizes plotted are: A'^ = 10 (up 
triangles), A^ = 12 (down triangles), A^^ = 14 (plus signs), and A^^ = 16 (right triangles). (Online version in colour.) 


as the average of the square overlap between many-body wave functions that evolved under perturbed Hamiltonians. 
While the former constitutes a specific experimental observable, the latter has only been assessed theoretically. Here, 
we showed the formal relation between both magnitudes, as far as their characteristic time scales and A^-dependence 
are concerned. 

By analyzing a short-time expansion of Midland Mmb we derived a precise the relation between their time scales. 
In this regime, the decay of is given by the average local second moment of the perturbation (fi/rs = and 

the decay of Mmb by N times the local scale {Nh/r^.). This relation hints a scaling law Mmb ~ that 
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accounts for the extensivity of Mmb- In such a case, the recovery of a many-spin state results from the recovery of 
each individual spin, much as if they were independent events. The numerical evaluation in a specific spin model shows 
that the exponent slightly diminishes with time, starting from the initial 7V/4. This means that the recovery of a single 
spin is positively correlated with the probability of recovery of its neighbors, and thus it improves the probability of 
the revival of the complete iV-spin state. A precise control of these correlations may hint an experimental access to 
the global autocorrelation, i.e. Mmb, just by measuring a single spin (local) autocorrelation M 14 . This would require 
an experimental protocol capable to encode a local excitation into a correlated many-spin state. 

In addition, we discussed a general dynamical picture beyond the very short-time regime. There, the decay of Mi 1 
results from the interplay between the time scale that characterizes the reversible interactions (T 2 ) and that of the 
perturbation (ts). This would ultimately lead to the experimentally observed T 3 , which was found to be roughly 
independent of rs but closely related to T 2 . The theoretical quest for the emergent T 3 time scale remains open and 
it may be out of the reach of current numerical approaches. Assessing a fair estimate analytically would require a 
detailed account of the higher order processes that dress the quadratic term in the perturbative expansion. 

Notice that our discussion lead us to identify T 3 , and hence the spin-spin interaction time T 2 , as the time scales 
characterizing the complexity or many-spin chaos. As such, they show up not only in the decay of Mi^i and Mmb, 
but also in the growth of Mx = Mi^ — Mmb- Indeed, in the field of AdS/CFT there is an increasing interest in 
characterizing the role of chaos in quantum dynamics [55l - l58l| . There, chaos manifests in the growth of four-body 
correlation functions, following an early suggestion by Larkin and Ovchinnikov (s^ . They employed semiclassical 
arguments to address disordered superconductors and probed that the square dispersion of momentum should grow 
exponentially in a time scale determined by the collisions with impurities, i.e. with the unperturbed Hamiltonian 
(in our physical picture, 72 ). Similarly, our average multispin correlation Mx would ultimately diverge within a 
time scale T^/N, i.e. independent of the perturbation. This is indeed a measure of the decoherence, and hence of 
irreversibility, induced by many-spin chaos. Of course, we do not have a precise characterization of this time scale 
or the specific mathematical dependence on time. Thus, this is a puzzling issue to explore in the field of many-body 
chaos. Besides the obvious relevance for statistical mechanics and experimental physics, this might also contribute to 
a possible pathway between quantum mechanics and gravity. 
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VI. APPENDIX 


It is worthy to mention that very short-time expansions in equations (1251) . (1261) and (1271) do not depend on the 


anisotropy a of the perturbation. In general, it can be proved that if 


^,si 


= 0 then: 
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(41) 


This is precisely the case of a perturbation E enclosing Anderson disorder and Ising interactions or interactions 
with a fluctuating field [^ . 


[1] Blume M, Hubbard J. 1970 Spin correlation functions at high temperatures. Phys. Rev. B 1, 3815-3830. 

[2] Forster D. 1975 Hydrodynamic fluctuations, broken symmetry, and correlation functions, Reading, Mass., WA Benjamin, 
Inc., Frontiers in Physics. Vol 47. 

[3] Basko DM, Aleiner IL, Altshuler BL. 2006 Metal insulator transition in a weakly interacting many-electron system with 
localized single-particle states. Ann. Phys. 321, 1126-1205. 












12 


[4] Zangara PR, Levstein PR, Pastawski HM. 2015 Role of energy uncertainties in ergodicity breaking induced by competing 
interactions and disorder, a dynamical assessment through the Loschmidt echo. Papers in Physics 7, 070012. 

[5] Wen XG. 2004 Quantum Field Theory of Many-body Systems from the Origin of Sound to an Origin of Light and Electrons, 
Oxford University Press Inc., New York. 

[6] Brewer RG, Hahn EL. 1984 Atomic Memory. Sci. Am. 251, 50-57. 

[7] Gorin T, Prosen T, Seligman TH, Znidaric M. 2006 Dynamics of Loschmidt echoes and fidelity decay. Phys. Rep. 435, 33 
- 156. 

[8] Jacquod P, Petitjean C. 2009 Decoherence, entanglement and irreversibility in quantum dynamical systems with few degrees 
of freedom. Adv. in Phys. 58, 67-196. 

[9] Goussev A, Jalabert R, Pastawski H, Wisniacki D. 2012 Loschmidt echo. Scholarpedia 7, 11687. 

[10] Hahn EL. 1950 Spin echoes. Phys. Rev. 80, 580. 

[11] Rhim WK, Kessemeier H. 1971 Transverse-magnetization recovery in the rotating frame. Phys. Rev. B 3, 3655-3661. 

[12] Rhim WK, Pines A, Waugh JS. 1971 Time-reversal experiments in dipolar-coupled spin systems. Phys. Rev. B 3, 684-696. 

[13] Zhang S, Meier BH, Ernst RR. 1992 Polarization echoes in NMR. Phys. Rev. Lett. 69, 2149-2151. 

[14] Levstein PR, Usaj G, Pastawski HM. 1998 Attenuation of polarization echoes in nuclear magnetic resonance: A study of 
the emergence of dynamical irreversibility in many-body quantum systems. J. Chem. Phys. 108, 2718-2724. 

[15] Pastawski HM, Levstein PR, Usaj G, Raya J, Hirschinger J. 2000 A nuclear magnetic resonance answer to the Boltzmann- 
Loschmidt controversy? Physica A 283, 166 - 170. 

[16] Usaj G, Pastawski HM, Levstein PR. 1998 Gaussian to exponential crossover in the attenuation of polarization echoes in 
NMR. Mol. Phys. 95, 1229-1236. 

[17] Flambaum VV, Izrailev FM. 2000 Excited eigenstates and strength functions for isolated systems of interacting particles. 
Phys. Rev. E 61, 2539-2542. 

[18] Flambaum VV, Izrailev FM. 2001 Entropy production and wave packet dynamics in the Fock space of closed chaotic 
many-body systems. Phys. Rev. E 64, 036220. 

[19] Peres A. 1984 Stability of quantum motion in chaotic and regular systems. Phys. Rev. A 30, 1610-1615. 

[20] Jalabert RA, Pastawski HM. 2001 Environment-independent decoherence rate in classically chaotic systems. Phys. Rev. 
Lett. 86 , 2490-2493. 

[21] Gucchietti FM, Pastawski HM, Jalabert RA. 2004 Universality of the Lyapunov regime for the Loschmidt echo. Phys. Rev. 
B 70, 035311. 

[22] Krojanski HG, Suter D. 2004 Scaling of decoherence in wide NMR quantum registers. Phys. Rev. Lett. 93, 090501. 

[23] Morgan SW, Oganesyan V, Boutis GS. 2012 Multispin correlations and pseudothermalization of the transient density 
matrix in solid-state NMR: Free induction decay and magic echo. Phys. Rev. B 86 , 214410. 

[24] Kaur G, Ajoy A, Gappellaro P. 2013 Decay of spin coherences in one-dimensional spin systems. New J. Phys. 15, 093035. 

[25] Sanchez GM, Acosta RH, Levstein PR, Pastawski HM, Ghattah AK. 2014 Glustering and decoherence of correlated spins 
under double quantum dynamics. Phys. Rev. A 90, 042122. 

[26] Alvarez GA, Suter D, Kaiser R. 2015 Localization-delocalization transition in the dynamics of dipolar-coupled nuclear 
spins. Science 349, no. 6250, 846-848 

[27] Zangara PR, Dente AD, Levstein PR, Pastawski HM. 2012 Loschmidt echo as a robust decoherence quantifier for many- 
body systems. Phys. Rev. A 86 , 012322. 

[28] Fine BV, Elsayed TA, Kropf GM, de Wijn AS. 2014 Absence of exponential sensitivity to small perturbations in noninte- 
grable systems of spins 1/2. Phys. Rev. E 89, 012923. 

[29] Zangara PR, Bendersky D, Pastawski HM. 2015 Proliferation of effective interactions: Decoherence-induced equilibration 
in a closed many-body system. Phys. Rev. A 91, 042112. 

[30] Pastawski HM, Levstein PR, Usaj G. 1995 Quantum dynamical echoes in the spin diffusion in mesoscopic systems. Phys. 
Rev. Lett. 75, 4310-4313. 

[31] Pastawski HM, Usaj G, Levstein PR. 1996 Quantum interference phenomena in the local polarization dynamics of meso¬ 
scopic systems: an NMR observation. Chem. Phys. Lett. 261, 329-334. 

[32] Madi ZL, Brutscher B, Schulte-Herbriiggen T, Briischweiler R, Ernst RR. 1997 Time-resolved observation of spin waves in 
a linear chain of nuclear spins. Chem. Phys. Lett. 268, 300 - 305. 

[33] Jacquod P, Silvestrov P, Beenakker C. 2001 Golden rule decay versus Lyapunov decay of the quantum Loschmidt echo. 
Phys. Rev. E 64, 055203. 

[34] Jacquod P, Adagideli I, Beenakker GWJ. 2002 Decay of the Loschmidt echo for quantum states with sub-Planck-scale 
structures. Phys. Rev. Lett. 89, 154103. 

[35] Prosen T. 2007 Ghaos and complexity of quantum motion. J. Phys. A-Math. Theor. 40, 7881. 

[36] Torres-Herrera EJ, Vyas M, Santos LF. 2014 General features of the relaxation dynamics of interacting quantum systems. 
New J. Phys. 16, 063010. 

[37] Torres-Herrera EJ, Santos LF. 2014 Quench dynamics of isolated many-body quantum systems. Phys. Rev. A 89, 043620. 

[38] Gappellaro P. 2014 Implementation of state transfer hamiltonians in spin chains with magnetic resonance techniques. In 
Quantum State Transfer and Network Engineering (ed. GM Nikolopoulos, I Jex), pp. 183-222. Springer Berlin Heidelberg. 

[39] Fukuhara T, Kantian A, Endres M, Gheneau M, Schaufi P, Hild S, Bellem D, Schollwock U, Giamarchi T, Gross G, Bloch 
I, Kuhr S. 2013 Quantum dynamics of a mobile spin impurity. Nat. Phys. 9, 235-241. 

[40] Zangara PR, Dente AD, lucci A, Levstein PR, Pastawski HM. 2013 Interaction-disorder competition in a spin system 
evaluated through the Loschmidt echo. Phys. Rev. B 88 , 195106. 



13 


[41] Serbyn M, Knap M, Gopalakrishnan S, Papic Z, Yao NY, Laumann CR, Abanin DA, Lukin MD, Dernier EA. 2014 
Interferometric probes of many-body localization. Phys. Rev. Lett. 113, 147204. 

[42] Levstein PR, Pastawski HM, Calvo R. 1991 Spin diffusion in low-dimensional copper-amino-acid complexes. J. Phys.: 
Condens. Matter 3, 1877. 

[43] Cucchietti FM, Dalvit DAR, Paz JP, Zurek WH. 2003 Decoherence and the Loschmidt echo. Phys. Rev. Lett. 91, 210403. 

[44] Bendersky D, Zangara PR, Pastawski HM. 2013 Fragility of superposition states evaluated by the Loschmidt echo. Phys. 
liev. A 88 , 032102. 

[45] Alvarez GA, Danieli EP, Levstein PR, Pastawski HM. 2008 Quantum parallelism as a tool for ensemble spin dynamics 
calculations. Phys. Rev. Lett. 101, 120503. 

[46] Elsayed TA, Fine BV. 2013 Regression relation for pure quantum states and its implications for efficient computing. Phys. 
Rev. Lett. 110, 070404. 

[47] Pineda G, Prosen T, Villasenor EV. 2014 Two dimensional kicked quantum Ising model: dynamical phase transitions. New 
J. Phys. 16, 123044. 

[48] Dominguez F, Gonzalez G, Segnorile H, Zamar R. 2015 Irreversible decoherence of dipole interacting nuclear spins coupled 
with a phonon bath. arXiv preprint:1505.01559 . 

[49] Anderson PW. 1978 Local moments and localized states. Rev. Mod. Phys. 50, 191-201. 

[50] Rodrigues DE, Pastawski HM, Weisz JF. 1986 Localization and phase coherence length in the Lloyd model. Phys. Rev. B 
34, 8545-8549. 

[51] Pastawski HM. 2007 Revisiting the Fermi golden rule: Quantum dynamical phase transition as a paradigm shift. Physica 
B 398, 278 - 286. 

[52] Lieb E, Schultz T, Mattis D. 1961 Two soluble models of an antiferromagnetic chain. Ann. Phys. 16, 407 - 466. 

[53] Rufeil-Fiori E, Sanchez GM, Oliva FY, Pastawski HM, Levstein PR. 2009 Effective one-body dynamics in multiple-quantum 
NMR experiments. Phys. Rev. A 79, 032324. 

[54] Dente AD, Bederian GS, Zangara PR, Pastawski HM. 2013 GPU accelerated Trotter-Suzuki solver for quantum spin 
dynamics. arXiv preprint:1305.0036 . 

[55] Shenker SH, Stanford D. 2014 Black holes and the butterfly effect. J. High Energy Phys. 3, 67. 

[56] Sekino Y, Susskind L. 2008 Fast scramblers. J. High Energy Phys. 10, 065. 

[57] Roberts DA, Stanford D, Susskind L. 2015 Localized shocks. J. High Energy Phys. 2015, 51. 

[58] Kitaev A. 2014 Hidden correlations in the Hawking radiation and thermal noise. Talk given at Fundamental Physics Prize 
Symposium, Nov. 10,. Stanford SITP seminars, Nov. 11 and Dec. 18. 

[59] Larkin AI, Ovchinnikov YN. 1969 Quasiclassical Method in the Theory of Superconductivity. J. Exp. Theor. Phys. 28, 
1200 . 

[60] Fernandez-Alcazar LJ, Pastawski HM. 2015 Decoherent time-dependent transport beyond the Landauer-Bfittiker formu¬ 
lation: A quantum-drift alternative to quantum jumps. Phys. Rev. A 91, 022117. 



